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ABSTRACT 

We study a cosmological scenario for gamma ray bursts (GRBs) where relativistic flows 
interact with dense radiation flelds. It is shown that this scenario is plausible in very 
dense stellar regions which are known to exist in collapsed cores of globular clusters 
or dense nuclei of galaxies. It yields a correct quantitative description of the temporal 
behavior of GRBs. Several other properties of GRBs are easily explained. 

Key words: globular clusters: general - gamma rays: bursts 



^1 INTRODUCTION 

Since their discovery, Gamma-Ray Bursts (GRBs) have re- 
^ mained an unsolved puzzle. Many of their features still pose 
I difficult problems. If originating at cosmological distances, 
as might be suggested by their remarkable isotropy (Meegan 
^^et al. 1992), then their production involves an enormous re- 
l/^ lease of energy in a short time. Their short durations imply 
0^very compact sources and enormous pair production, which 
makes them optically thick to their own 7-ray emission. 
H ence, the original 7-rays are concealed by an opaque enve- 
I lope, which reduces the original 7-rays to much softer ther- 
O mal photons (Cavallo & Rees 1978, Paczynski 1986, Good- 
-j-J man 1986, and Goodman, Dar & Nussinov 1987). They are 
^ emitted when eventually the fireball expands enough for the 
optical depth to decrease below unity. Although the expan- 
sion does not necessarily increase the duration because of its 
Lorentz beaming of the emitted radiation, it cannot explain 
the non-thermal spectrum of GRBs. Furthermore, baryonic 
contamination greatly affects its dynamics. Even a baryon 
load as small as Mbaryon ~ 10~^Mq can convert most of 
the energy stored in the radiation to kinetic energy of the 
baryons (Paczynski 1990, Shemi & Piran 1990). 

It turns out, that both the non-thermal spectrum and 
the effects of the baryon overload can still agree with fire- 
ball scenarios if the kinetic energy of the baryons is con- 
verted back into radiation. Several such mechanisms exist. 
For instance, Rees & Meszaros (1992) suggest the decelera- 
tion of the relativistic wind by the ambient interstellar mat- 
ter, which is compressed in relativistic shock waves and then 
emits 7-rays while cooling. An alternative mechanism for 
reconversion of the kinetic energy into 7-rays was proposed 
by Shemi (1993a) where the interaction of the relativistic 
wind is with a local thermal radiation field. The relativistic 
electrons up scatter the ambient photons, and shift them 
towards the direction of flow, modifying their spectrum into 



a power law. The rather large scale of ^ 10'^ light-seconds 
is still transformed into a short burst because of the rela- 
tivistic beaming. Furthermore, if the wind is not in a form 
of a fireball, but of a jet, then the large spectrum of dura- 
tions can easily be explained as a geometrical effect (Shaviv 
& Dar 1995). In his work, Shemi (1994) has shown that 
Compton drag of relativistic debris of cosmic fireballs can 
explain several generic problems in optically thick GRBs - 
their compactness, non-thermal spectrum, effects of bary- 
onic loads, and time scales. It has though one very major 
drawback - it needs very dense radiation fields, those occur- 
ring only in dense cores of globular cluster and in galactic 
nuclei. 

We first consider the scenario proposed by Shemi 
(1993a, 1994) of a thin fireball shell (as is shown to form by 
Shemi 1993b, Piran, Shemi & Narayan 1993 and Meszaros, 
Laguna & Rees 1993) scattering photons of the radiation 
field in a dense stellar cluster. We then show that in addi- 
tion to explaining various properties of GRBs, it yields time 
profiles which agree remarkably well with the experimental 
observations. 



2 ENERGETICS OF RADIATIONALLY 
DECELERATED FIREBALLS 

How plausible is it to find high photon density regions? One 
has to go only a few kpc's to reach the collapsed cores of 
several globular clusters (about 20% of the globular clusters 
are core collapsed, see e.g. Chernoff & Djorgovski 1989) or 
the centre of our galaxy. In several globular clusters which 
had their core collapsed, one can find within a sphere of 
radius 0.1-pc a large fraction of the cluster's mass, giving 
very high stellar densities (exceeding even lO^A^0pc"^ e.g., 
Phinney & Sigurdsson 1991). Unfortunately, the centre of 
our galaxy is blocked by gas and dust, but in the galaxy 
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of Andromeda, where direct visible measurements can be 
made, there exists a region of 5-pc by 10-pc which contains 
a total mass of 10'' to 10* corresponding to a stellar 

density of about 2 x 10 Mqpc~ . As we shall see, these 
types of stellar densities can give rise to an ample photon 
field, capable of transforming a great fraction of the kinetic 
energy of fireballs into 7-rays. 

It is interesting to note that the dynamics of collapsed 
cores of globular clusters must give rise to hard binaries 
which stop the unstable collapse (e.g. Heggie 1980). These 
abundant binaries are the probable cause for the relatively 
very large number of 'interesting' objects such as low mass 
X-ray binaries and binary neutron stars. Among the four 
known systems of neutron star binaries, PSR 2127+110 is 
in the globular cluster M15 (Prince et al. 1991), in contrast 
with the total mass in globulars which is only 10"'^ of the 
total galactic mass. It is thus very plausible that these dense 
clusters hold within them many GRB progenitors, much 
more than can be estimated from their mass fraction ^ . 

If we look at such a dense cluster having a stellar pop- 
ulation of within a region of radius R = 0.1 rfo.ipc (do.i 
measures the radius in units of 0.1-pc), and an average stel- 
lar luminosity ILq with photons of energy eo-eV, then the 
photon density at the centre is: 



o n 1 •^^*.5 -3 

2.0 X 10 — — 7 cm , 



(1) 



with = 10^ ^★.B. If we assume such a photon density, and 
that the only means of interaction is by Compton scatter- 
ing, then we can calculate the amount of energy lost by an 
electron while traversing a path of length R (i.e. the average 
distance until it leaves the cluster). 

The electron's energy is: Ee = frUeC^ ■ On the average, 
an electron loses E-y = jf^^o of its energy per collision. The 
mean free path is A = l/crTn-^ where ctt = 0.66 x 10~^*cm^ 
is the Thomson cross section, and the average number of in- 
teractions along R is R/X = Rutu.^. Therefore, the electron 
will lose an average energy of eoRaTti-f before leaving 
the cluster. Its fractional energy loss, ^, is given by: 

4(7Tra7-R7eo 
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(2) 



with meff being the effective mass associated with each elec- 
tron. For a baryonic fireball, meff = me -|- nip K, nip, but 
if the fireball is still rich in pairs, the effective mass can be 
lower, reaching nie for Wpairs >> 1000 Wp. 

Since the duration of the GRB is given by to ~ 
7~^(i?/2c), one obtains: 



c 



4.0 X 10"*'£7V^,5V^^^(i^y^/tt"\ 



(3) 



with fi = meff/mp. For the Andromeda galaxy where iV^_5 Si 
1000, (i^y^ f» 0.1 and t~^''^ = 0.3 or 3 for 10 sec or 0.1 sec 
bursts, eq. (3) for /t = 1 yields C = 0.1 x 10"^ or 1 x 10"^ 
respectively. If the observed 7-ray energy is £50 x 10^"- 
erg then the initial kinetic energy of the fireball, assum- 
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Figure 1. A schematic drawing illustrating the proposed sce- 
nario. A fireball inside a dense stellar region is formed. While 
expanding it scatters the cluster's radiation field photons. Most 
of the observed scattered photons come from the region near the 
stars, were the light density is the highest. 



respectively. For collapsed cored globular clusters (such as 
M15) we have smaller values: 7V^_5 Si 5 and dg K, 1 
yielding: C, = 0.5 x 10"'' or 5 x 10"'' for /t = 1, giving 
-Ei£in = 1000i?5o or 100i?5o Mqc^, respectively. These num- 
bers are obviously too large, rendering this model impossible 
unless we either utilize beaming, decrease fi or increase a. If 
we succeed to form a pair-rich fireball, fi can reach 2000, in 
which case the required initial kinetic energy can range from 
-Eitin = 2.5 X 10"*_£:5o Mqc^ to 0.5_£:5o Mqc^ , depending on 
the parameters, without any beaming. 

Having a higher cross-section a for scattering photons 
off the fireball is a very realistic possibility for NS-NS merger 
scenarios. It is well known that the crust of neutron stars is 
composed of mainly ^^Fe or heavier nuclei. A fireball com- 
posed of such nuclei will have a very large opacity at a tem- 
perature of a few eV ^. For such nuclei, this opacity (per 
unit mass in the fireball) formed by the large cross-section 
for ionization by the ambient photons that in the fireball's 
rest frame have keV energies, can be greater than 10^ times 
the opacity resulting from Compton Scattering. The ionized 
levels are quickly filled, emitting back photons that in the 
observers frame appear as ~ 100-keV photons. Since the dif- 
ferent lines are so numerous, the smearing resulting from the 
different Doppler shifts of the different shells in the fireball 
will smear out any possible lines. Here, the required initial 



kinetic energy can range from -Eitin 



5 X 10"^ £50 M, 



to O.Ol^Bo Mqc\ The beaming geometry will reduce the 
energies even more. 



ing no beaming, is i?i£in 



50_£:5o Mqc^ or 5_£:5o M, 



0' 



^ For Instance, the conservative estimate by Phinney (1991) 
of the neutron star merger rate in globular clusters is 3 X 
10~^/i'^Mpc~'^yr~^. This rate is smaller by less than an order 
of magnitude than the rate required to meet the CGRO-BATSE 
event rate 



^ The fireball's temperature can be found by calculating the tem- 
perature and radius at which the photon opacity has decreased to 
unity. Such a calculation (with equations that can be found, for 
example, in Meszaros, Laguna Rees 1993) yields a temperature 
of a few eV. At these temperatures and densities, all atomic levels 
with energies less than ~ 200-eV will be ionized. 
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3 THE TEMPORAL PROFILE OF A 

FIREBALL INTERACTING WITH DENSE 
INTERSTELLAR RADIATION 

We shall now proceed to calculate the time profile of a GRB 
resulting from a fireball shell that substands a large opening 
angle and expands in a dense stellar field. If we look at Fig 1, 
we see that two regions can emit most of the observed GRB 
energy. The first is the region along the observer's line of 
sight, while the second is the region near a star. The former 
is a region of angular size x l/~f^ with a boosting factor of 
X 7, while the second is a region of angular size 8^ with a 
boosting factor of: 



1 



1 



7612 + I/7 7612 



(4) 



leading to the same result. However, these factors must be 
weighted with the appropriate distance attenuation, which 
is larger for the former region. Hence, the received fluence 
will come mainly from a region near the star. We shall now 
estimate its time profile. The time dependent intensity pro- 
file formed by just one star is given approximately by: 
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If J'min which is formed from the finite width of the shell 
is small, then between a;min which is of the order of one over 
the duration of the burst and a;max, there will be a large 
region where /i will have a 1/a; behavior. The total power 
spectrum in this region is given by: 



Ptotal(t^) 




E 



(8) 



Since one has to sum different Fourier transforms with ran- 
dom phases one obtains the second equality, which simply 
states that the power spectrum of a sum of uncorrelated 
Fourier transforms is just the sum of the different individual 
power spectra. Note that the total power spectrum in the 
middle region has a behavior, which as we shall see 

later, is the power spectrum behavior of real GRBs. This 
l/uP behavior will break down for small enough a;'s, smaller 
than a few times the reciprocal of the duration, and also for 
large enough a;'s. If the shortest features in GRBs are of 
10~'^s, and originate from the finite width of the shell, than 
the upper frequency at which the law breaks is ~ 1000 Hz. 
A larger shell size will result with a lower upper frequency. 



where (7c(^)rfQ is the fraction of radiation scattered into 
the solid angle rfQ at an angle of 9 from the observer's 
line of sight. This cross-section approximately behaves as 
oc 1/(7^^ -|- 1/7). rmin is a lower cutoff that is formed by the 
finite width of the shell, and it is assumed to be much less 
than fmax. ti is the time when the shell crosses the star, and 
Io,i is the intrinsic luminosity of the star. Note that this re- 
sult is a very rough approximation. Eq. (17) which assumes 
that different parts of the fireball scatter with the same in- 
cident angle 9, essentially approximates the fireball with a 
round plane. This assumption is justified since most of the 
scattered radiation comes from a small region of the fireball 
near the star, which does not see the curvature of the fire- 
ball. In addition, it does not take into account the fact that 
the scattered intensity depends on both the initial spectrum 
of the photons and the scattering angle. These effects are 
too complex to be taken into account analytically. Together 
with the finite width of the fireball they probably tend to 
deform the resulting time profile. 

The Fourier transform of eq. (17) is: 



/i(o 



^ / I.{t)e"^'dt 

ZTT / 

^ —CO 



dt 



4a; 



(6) 



It has three distinctive regions of behavior bounded by 
o^min = f/^max and a;max = v/ Train- The bchavior in each 
region is given by: 



4 NUMERICAL SIMULATIONS AND 

COMPARISON WITH OBSERVATIONS 

In addition to the analytical calculations, we have performed 
several numerical simulations of fireball expansion in dense 
stellar regions. Although the analytical calculations yield a 
temporal behavior similar to that observed for GRBs, they 
contained several approximations which need verification. In 
addition, the numerical simulations give a general picture of 
the emerging time profiles, which is usually not evident from 
just the power law behavior. 

In Fig. 2 we present a typical temporal behavior cal- 
culated for a Gaussian star cluster with = 5 x 10^ and 
with a radius of do.i = 1, fmin = 10~*rfo.i and 7 = 100. 
The temporal profile is normalized to the burst's length of 
1/7^. The power spectrum of the temporal profile in Fig. 2 
is depicted in Fig. 3. Fig. 4 and 5 are the same as Fig. 2 and 
3, but with addition of simulated noise. In both power spec- 
tra, the law is evident. In the second case, the power 
slope can be seen down to frequencies which are lower than 
the frequency where the signal is below the noise. Figures 
6 and 7 show the temporal profile and the power spectrum 
of an arbitrarily selected real GRB taken from the BATSE 
GRB catalogue (Fishman et al. 1994). The simulated pro- 
file and the real profile seem to have the same 'qualitative' 
and 'quantitative' behaviors. That is the case also for other 
GRBs since the simulated burst's spectral power law be- 
havior is the same as for real GRBs (Shaviv & Dar 1995) - 
namely a behavior. 

Fig. 8 and 9 show the temporal and spectral result of 
a Gaussian star cluster with = 10^ and with a radius of 
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Figure 2. The temporal profile of a simulated GRB from a Gaus- 
sian star cluster with A^^ = 5 X 10^ and with a radius of unity. 
In addition rmin = 10~^ and 7 = 100. The temporal profile is 
normalized to the burst's length of 1/7^. 
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Figure 3. The power spectrum of the time profile shown in Fig. 
2. 



do.i = 1, but with 7 = 1000. Since the fireball seems more 
squeezed to an outside observer, less stars will pass through 
it, giving less peaks inside. Note that the simple profile of 
eq. 17 was used. More realistic assumption on the scattering 
process and the finite width of the fireball will give more 
complex asymmetric peaks, as can be seen in the real GRB 
in Fig. 10 (with its Fourier transform in Fig. 11). Actual 
intensities will depend on the spectral band of the detec- 
tor and the spectrum of both the scattering electrons and 
photons, which together with the scattering angle change 
the outcoming energy of the photons. This effect is proba- 
bly negligible, since it is a slowly varying function of time, 
giving a Fourier transform which is 'localized' in frequency 
space, and therefore does not change the law when 

convolved with the previous result. 

The model can explain several additional observations. 
For instance, there have been several reports of after burst 
pulses of very energetic photons (e.g. Mukherjee 1994, Hur- 
ley et al. 1994). These can come from an X-ray source far 
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Figure 4. Same as Fig. 2, but with a small white noise component 
added to simulate measured profiles. 
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Figure 5. The power spectrum of the time profile in Fig. 4. Two 
lines are drawn, one shows an uj~^ fit, while the other shows the 
noise level. 
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Figure 6. The temporal profile of GRB920110 (trigger #1288 
of CGRO-BATSE). 
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Figure 7. The power spectrum of the time profile shown in Fig. 
6. 
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Figure 9. The power spectrum of the time profile shown in Fig. 
8. Two straight lines are drawn, one is the noise level, while the 
other has a slope of uj~^ . 
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Figure 8. The temporal profile of a simulated GRB from a Gaus- 
sian star cluster with A^^ = 10^ and with a unit radius. 7 = 1000. 
The temporal profile is normalized to the burst's duration of 7"^ . 
A small white noise component was added. 



from the observer's line of sight, where the Lorentz beaming 
is too weak to boost the optical radiation into 7-rays, but 
can boost the X-rays into hard 7-rays . Since the source is 
far from the line of sight, the time for the fireball to reach 
this X-ray source, as measured by the observer, will be much 
longer. An additional possibility is the annihilation in flight 
of positrons left in the fireball. The number density of elec- 
trons in the interstellar matter in the galaxy, and certainly 
in globular clusters, is much lower than the number density 
of the optical photons in the centre of the cluster. As a con- 
sequence, we can expect that only a small fraction of the 
positrons would annihilate in flight while interacting with 
electrons of the interstellar medium. These events could also 
occur out of the dense stellar regions, resulting with delayed 
hard photons. 



X-ray sources are known to be abundant in globular clusters, 
see e.g. Katz (1975) and Lewin (1980). 




Figure 10. The temporal profile of GRB911031 (trigger #973 of 
CGRO-BATSE). 



Ford et al. (1995) have reported that there is a general 
softening of the GRBs with time. Occasionally though, the 
bursts harden with the increase of the intensity. This also 
coincides qualitatively with the model. As the wind is de- 
celerated, the lower Lorentz 7 factor yields on the average a 
softer spectrum. But, when passing nearby a star, the hard- 
ness depends upon the relative incident angle. Occasionally 
a brighter peak and a harder spectrum will be produced 
while passing a star. 



5 DISCUSSION 

At first sight, the interaction of a fireball with dense radia- 
tion fields seems highly improbable as the source of GRBs. 
However, we have shown that fireball debris can lose enough 
of their energy and generate a GRB if occurring in dense 
stellar fields. Such fields exist in dense galactic nuclei and in 
the cores of collapsed globular clusters. As was shown, this is 
possible if we can produce a fireball containing many e"'"e~ 
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Figure 11. The power spectrum of the time profile shown in Fig. 
10. The two straight lines are again the noise level and an uj~^ 
fit. 

pairs left after freeze out (as compared witli tlie number of 
baryons) which do not annihilate. This enables the photons 
from the radiation field to slow down the massive debris, 
composed mainly e"'"e~ pairs. An additional possibility ex- 
ists, if the fireball's opacity to the ambient photons is large. 
This results naturally in fireballs that are composed of high 
Z material which can originate, for example, from the crust 
of neutron stars. A problem would occur though, if the in- 
terstellar medium is not rarefied. A medium containing too 
much gas (wism ^ 1 cm~'^) can possibly stop the fireball 
without releasing the energy by boosting the radiation. 

The structure of the dense radiation fields was taken to 
be that formed from a random distribution of stars. This 
was used to calculate the power spectrum of the temporal 
profile, and gave a power law spectrum with a power index 
of —2 - exactly the same power spectrum of the observed 
GRBs. It is interesting to note that even if the field contains 
only a single star, it will produce the same power spectrum; 
i.e., the power spectrum is not due to the statistical average 
of many stars, but it is a feature of the individual peeks. 

It is clear that the model is still in its infancy, but it 
already yields important verified predictions. Yet, there are 
still various GRBs features that have to be explained by the 
model. One such feature is the spectra of GRBs, which in 
the model, depends on the average spectrum of the scattered 
radiation, and the distribution of Lerentz boosting factors 
in the fireball. 
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